Stability mechanism of cuboctahedral clusters in U02 +x : First-principles calculations 
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The stability mechanism of cuboctahedral clusters in nonstoichiometric uranium dioxide is in- 
vestigated by first-principles LSDA+U method. Calculations reveal that the structural stability 
is inherited from UeOi2 molecular cluster whereas the energy gain through occupying its center 
with an additional oxygen makes the cluster win out by competition with point oxygen intersti- 
tials. Local displacement of the center oxygen along (111) direction also leads the cluster 8- folded 
degeneracy and increases relatively the concentration at finite temperatures. But totally, elevation 
of temperature, i.e., the effect of entropy, favors point interstitial over cuboctahedral clusters. 
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Uranium dioxide adopts the simple fluorite (CaF 2 ) type of crystal structure with a space group Fm3m. But its 
self-defects, as in most anion excess fiuorites, exhibit rather complex behavior: experimentally oxygen interstitials 
do not occupy the largest cation octahedral hole (-,-,-) but form low symmetric clusters composed of oxygen 
vacancies insteadi^ The exact geometry of these clusters, however, is unclear. Experimentalists have proposed 
several structural models to explain the measured neutron diffraction patterns^ of which the cuboctahedral cluster 
(COT) appearing in JJ^Og^^ and U3O7 7,8 is the most clearly described. But even with this structure, ambiguity 
remains about whether the center is occupied by an additional oxygen (COT-o) or not (COT-v), 7 and if the center 
atom really displaced off-center along (111) direction in the case that it was occupied. 51 - Theoretical analysis have 
done little help so far since most of work confined to point oxygen interstitial (0$) and vacancy (O v ) where the former 
always sites at the octahedral hole and thus failed to explain the experimental phenomena .— i 1 1 1 1 Inspired by its close 
relationship with Willis type clusters, COT was suggested also should present in U02+^. 4,9 A fully understanding 
of its property and stability mechanism thus becomes important not only for a general description of fluorite-related 
clusters 12 but also for nuclear applications, for example the safety disposal of used nuclear fuel where the formation 
of U4O9/U3O7 is a key process for the development of U3O8 phase which can lead to splitting of the sheath J£ 

Historically, COT was denoted by M6X36 (or M6X37 if an additional anion occupies the center) £^ It is a little mis- 
leading since the actual defect composed o f only 12 interstitials [forming the cuboctahedral geometry, see Fig l(a)| and 



8 vacancies [forming the small cube in Fig l(a)| . Different from Willis type clusters whose stability can be interpreted 



in a similar concept of split-interstitial defect where several atoms share the common lattice sites,- COT is of more 
regular and with higher symmetry [point group Oh(m3m)], and poises as the special one. In fact, our calculations 
that will be reported here revealed it is actually an U6O12 molecular cluster incorporated in bulk fluorite UO2 by 
sharing the uranium atoms with the cation sublattice face centers [Fig.l(a-b)] after removed the eight corresponding 
lattice oxygens from the matrix. The thus inserted 12 oxygens presented as Willis O interstitials that displaced along 
(110) directions in the picture of fluorite structure. 

In calculations, COT cluster was modeled by embedding it into a cubic supercell of fluorite UO2 with otherwise 96 
atoms (U32O64). This configuration has large enough cell size with a deviation composition x = | for COT-v and ^ 
for COT-o, compared with experimental x ~ 0.21 for U4O9 6 and 0.35 for U3O7, 7 respectively. U6O12 molecular cluster 
was modeled by put into a vacuum cubic box with a lateral length of 11 A, a sufficient distance for current purpose. 
Total energies were calculated with plane wave method based on density functional theory (DFT), 14 with generalized 
gradient approximation (GGA, for U6O12), local density approximation with Hubbard correction (LSDA+U ^ 15 i 16 for 
COTs) and projector-augmented wave (PAW) pseudopotentials i 17 i 18 All structures have been fully relaxed to get all 
forces and stress less than 0.01 eV/A. Other computational parameters such as energy cutoff and sampling k-points 
are the same as those in Ref.[9| which focused on point defects behaviors. 

Geometry and energetics. — The uraniums in COTs that embedded in UO2 are found always bond to interstitial 
oxygens firstly with a shorter bond length than to the nearest neighbor (NN) lattice oxygens (2.2 vs 2.4 A). Analysis 
of electron ic den sity also shows the weak covalent bonds that forming an U6O12 cluster are always prior to other 
bonds [Fig |l(d)] . Local distortions have not changed the picture and preserve much of the tightly connected feature 
of vacuum U6O12 cluster, which consists of three mutually perpendicular U-0 rings that in turn determined b y the 
two radius from the cluster center to the uranium (ru) and oxygen (ro) atoms completely, as shown in Fig l(c)| 
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FIG. 1: (Color online), (a) Cub octahedral cluster (COT-v) incorporated in a fluorite cell, where the small inner cube indicates 
the (removed) oxygen cage; (b) U6O12 molecular cluster; (c) One of the three mutually perpendicular U-0 rings in U6O12 
cluster, and (d) its charge density. 



TABLE I: First principles results for structural and energetic properties of oxygen defects in uranium dioxide: x is the deviation 
from the stoichiometric composition, AV the defect induced volume change that averaged to per fluorite cubic cell; ru and 
ro the structural parameters of cuboctahedral cluster as indicated in Fig |l(c)| L, UOU and OUO are the averaged nearest 
neighbor bond length between U and O atoms, and the corresponding angles in cuboctahedral cluster, respectively; E/, E e / 
and E p f are the overall defect formation energy, formation energy per excess oxygen and per oxygen Frenkel pair, respectively. 
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a Experiment of /3-U4O9 at 503 K reported in Ref.[6] 
b Isolated point defects, Ref.Q 



Vacuum U6O12 is perfect and without any deformation freedom such as variation in U-0 bond length and related 
angles. All real but low vibrational frequencies (not shown here) indicate the structure is locally stable but flexible 
for distortion. A cohesive energy of 21.2 eV per UO2 molecule is comparable with the bulk material [22.3(23.9) eV 
of experiment (GG A)] ^ These features and the geometry make U6O12 can be incorporated naturally with fluorite 
crystal (or generally, FCC lattice) and forms COT clusters without disturb the host structure severely. 

Locally, COTs repulse the NN lattice oxygens outwards slightly, but no evident distortion on cations was observed. 
As listed in table H the overall volume is contracted but the cube that contains the COT is expanded greatly, with a 
lateral length of 2r\j. Also, embedding U6O12 into UO2 not only swells the cube (with an increased ru and ro) but also 
leads to other local distortions. Usually the distorted U-0 ring is not on the same plane any longer, and has additional 
freedoms in U-0 bond length (L) and related angles (UOU and OUO). Their averaged values are listed in table U 
by compared with vacuum U6O12 and experimental estimates of COT measured at 503 K on /3-U4O9. 6 Obviously, 
COT-o agrees with the experimental one much better than COT-v in geometry. With an additional oxygen occupied 
the center, COT-o decreases the values of ru and UOU while lifts ro and OUO significantly with respect to COT-v, 
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FIG. 2: (Color online). Defect concentrations of point oxygen interstitial, uranium vacancy and COT-o cluster, respectively. 
Dashed lines indicate the corresponding results where the center oxygen in COT-o has no off-center displacement and thus with 
= 1- 



identifying it as the one that appeared in experiments. 

In contrast to the symmetry anticipation, 5 the center oxygen in COT-o does not site at the real center. It displaces 
along (111) direction of fiuorite structure with a distance about 0.43 A, compared well with Cooper and Willis' 
estimate of 0.64 A, 6 and stands as an O interstitial. In vacuum U6O12 cluster, such kind of off-center displacement is 
forbidden. The small r\j ensures the energy minimum always at the cluster center. But the r\j of COTs are enlarged 
by bulk UO2 matrix which in turn shifts the energy minimum off center. By displacing the center oxygen along 
(111) direction, the three NN uraniums within the corresponding section of the COT shell are pulled inwards ~ 0.2 A 
(while the three oxygens are pushed out slightly ~ 0.1 A) and reduced the U-0 bond length from 2.76 to 2.44 A. 
The Oh{m3m) symmetry is also broken to C3 V (3m). In contrast, the symmetry broken in COT-v is mainly due to 
Jahn- Teller distortion, where one uranium atom out relaxed but another five shift inwards, results in a C4 V (4rara) 
symmetry. 

The large cohesive energy of U6O12 leads to a deep formation energy of COT clusters, as indicated in table H where 
the energetic information of isolated point defects (0$ and O v ) are also given for comparison. By compensating excess 
oxygens with point vacancies, we find the formation energy per Frenkel pair in COT is just one third of the isolated 
case. However, the energy gain for each excess oxygen exhibits different behavior for COT-v and COT-o. With the 
contribution of the center oxygen, the latter has a lower E e / than the point interstitial but the former is weighed 
down by the electronic density cavity presented in the cluster center which costs the energy significantly. 

Defect concentrations. — With regard to the concentration of COT at finite temperatures, an intuitive picture 
is that it should favor moderate temperatures since otherwise there have not sufficient vacancies to facilitate the 
formation of the cluster. However, as mentioned above, the vacancies in COT actually are not from point defects but 
inherited integrally from the U6O12 molecular cluster, this naive picture thus becomes invalid. To calculate the defect 
concentrations at finite temperatures properly, we employed here the independent clusters approximation (ICA) (a 
generalization of the point defect model) ^ 2Q i 21 in which all involved clusters are assumed to be thermodynamically 
independent and obey Boltzmann distribution. In closed regime where no particle-exchange with the exterior occurs, 
the concentration pi of cluster i that has an internal degeneracy g^, m excess oxygens and a formation energy is 
given by 

f-Ei - mx E° v \ 
Pi\y r=9 i e X p f f — f j . (1) 

Here rii point oxygen vacancies (O v ) have been introduced as compensations. Therefore we have a system contains 
point oxygen(uranium) defects (treated as intrinsic Frenkel pairs), COT-v and COT-o clusters that compete to 
each other. The oxygen and uranium subsystems are coupled up via the isolated point Schottky defects £ Among 
these defects, only COT-o has an internal freedom with 8-folded degeneracy (g = 8) arising from the (111) direction 
displacement of the center oxygen and all others have g = 1. COT-v has 4 and COT-o has 5 excess oxygen, respectively. 
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FIG. 3: (Color online). Defect concentrations at different temperatures: solid (dotted) line — oxygen interstitial(vacancy) arising 
from COTs, dashed line — point oxygen interstitial, and dash-dotted line — point uranium vacancy. All other components are 
negligible. 



At a composition x deviated from the stoichiometry, Eqs. (pQ) are under a constraint of 



where the quantities in brackets denote point defect concentrations and i runs over COT-v and COT-o clusters, with 
a coefficient equals to the corresponding Ui because each octahedral hole of the cation sublattice defines not only a 
point interstitial site, but also a COT-v(o) cluster. 

Using the calculated first-principles formation energies, we get the defect concentrations as a function of temperature 
and composition by solving Eqs. (pQ) and ([2|). The hypostoichiometric regime is always dominated by O v and thus 
trivial. Interesting competition between defects appears on the other section of composition with x > 0. Solid lines 
in Fig. [2] show the equilibrium concentrations of point oxygen interstitial O;, uranium vacancy XJ V and COT-o cluster 
at a temperature of 1500 K around this regime. All other defects have a concentration of ten more orders smaller 
and not shown here. 0$ is predominant at low x with a rapid increment of its concentration, which flattens out 
gradually at high compositions. By contrast COT-o approaches a linear dependence on composition after x > 0.025 
and dominates the regime of x > 0.1. The overall concentration of XJ V is one order smaller than Oi and has dismissed 
its influence on material properties. 

The predominance of COT-o over COT-v is due to the energy gain of the center oxygen but not the entropy 
contribution of the 8-folded state. The dashed lines in Fig. [2] give the corresponding concentrations by treated COT-o 
with g = 1. Removal of the internal entropy contribution do reduce the competitiveness of COT-o and increase the 
concentration of Oi greatly, as well as that of XJ V , but has not changed the picture qualitatively — the concentration 
of COT-v is still tens more orders smaller. Conversely, if rescale the formation energy of COT-o so that it has the 
same E e f as COT-v, then they will have comparable concentrations, but of ten orders smaller than that of Oi. 

Figure [3] shows the variation of defects competition along temperature from 500 to 1500 K. The solid(dotted) line 
indicates the oxygen interstitial(vacancy) concentration arising from COT clusters: each COT-o(v) contributes 13(12) 
oxygen interstitials and 8 vacancies. It is obvious that increase temperature, i.e., the entropy effect, favors point defects 
over COT clusters greatly. The predominant range of Oi has increased 3 orders by elevate the temperature from 500 
to 1500 K. All point defect concentrations are enhanced along this process except those from COT clusters, which are 
reduced by temperature. This is because the probability to form a COT-o cluster from point defects is proportion to 
[7o] 13 [Vb] 8 but each COT-o in conversely contributes only 13(8) interstitials (vacancies), showing point defect is more 
disordered and with larger entropy. 

Summary and Discussion. — Figure [3] demonstrates that at 500 K COT-o is the exclusive defect cluster, support the 
empirical assumption that /S-XJ^Oq contains only this kind of cluster & From the temperature dependence of defect 
concentrations, we can be sure that the lower temperature a phase also should contain COT-o exclusively. On the 
other hand, the current interpretation of the neutron diffraction pattern in U3O7 is questionable, which employed 
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COT-v instead of COT-o cluster.— The former has a negligible concentration of ten more orders smaller and thus 
invalidates the analysis definitely. The ordering of the clusters that distributed in these phases^ however, seems 
should be driven by long-ranged strain energy rather than by chemical interactions. By the volume change induced 
by defects listed in Table U and the fact that COT itself expands the occupied fluorite cube seriously, there is a strong 
deformation field around each COT cluster which repulses other COTs away. The stress magnitude can be estimated 
from the bulk modulus of UO2 as ~ 2 GPa, a high enough value and any off-balance happened on the boundaries of 
deformed domains will lead to cracks. That is why U4O9/U3O7 film cannot protect UO2 pellet from being oxidized 
effectively^ Such cracks are also believed as the onset of high burn-up structures in nuclear fuels where uraniums are 
highly consumed and deteriorates the fuel quality severely. 22 
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